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We experimentally demonstrate the effect of anomalous breakdown of the effective medium ap¬ 
proximation in all-dielectric deeply subwavelength thickness (d ~ A/160 — A/30) multilayers, as 
recently predicted theoretically [H.H. Sheinfux et ak, Phys. Rev. Lett. 113, 243901 (2014)]. Mul¬ 
tilayer stacks are composed of alternating alumina and titania layers fabricated using atomic layer 
deposition. For light incident on such multilayers at angles near the total internal rehection we 
observe pronounced differences in the rehectance spectra for structures with 10-nm versus 20-nm 
thick layers, as well as for structures with different layers ordering, contrary to the predictions of 
the effective medium approximation. The rehectance diherence can reach values up to 0.5, owing to 
the chosen geometrical conhguration with an additional resonator layer employed for the enhance¬ 
ment of the ehect. Our results are important for the development of new high-precision multilayer 
ehipsometry methods and schemes, as well as in a broad range of sensing applications. 

PACS numbers: 78.67.Pt, 78.20.-e, 42.25.Gy, 81.15.Gh 


Photonic multilayers are one of the most widely stud¬ 
ied systems in the broader topic of optics of inhomoge¬ 
neous media M. Most optical effects of such multilay¬ 
ers arise from interference effects underlying the photonic 
band gap phenomena @ , and therefore are traditionally 
associated with multilayers, where layer thickness d is 
comparable to the wavelength of light A. For example, 
the well-known Bragg mirror comprising a stack of al¬ 
ternating low- and high-index dielectric layers, exhibits 
maximum reflectance if the optical thickness of each layer 
is close to A/4 0. 

From this point of view, the case of multilayers with 
much thinner layers with thicknesses d X was tra¬ 
ditionally regarded as nearly trivial. Indeed, in such a 
“deeply subwavelength” structure the field variation in¬ 
side a single layer should be very small, leading to neg¬ 
ligibly weak interference effects. Therefore, one used to 
assume that a light wave interacts with the structure as 
a whole rather than with its individual layers. The struc¬ 
ture can thus be treated as a piece of homogeneous uni¬ 
axial material characterized by effective parameters 0- 
The applicability of this homogenization approach to all¬ 
dielectric multilayers with ultrathin layers d X has al¬ 
ways been undoubted, in much the same way as ordinary 
materials are treated as homogeneous media despite hav¬ 
ing atomic, molecular, or any other intrinsic structure. 

However, a recent paper by Sheinfux et al. 0 showed 
theoretically that this commonly believed assumption 
may fail in certain circumstances. Namely, when light is 
incident on a multilayer at an angle close to that of the 
total internal reflection (TIR), the actual multilayer and 
its effective-medium approximation (EMA) model can 


have significantly different transmission spectra despite 
the layer thicknesses smaller than A/50. Moreover, it was 
shown that the spectra become sensitive to variations of 
d on the scale of one nanometer (i.e., < A/500), as well 
as on the layer ordering (i.e., the multilayer HLHL...HL 
versus LHLH...LH, where H and L stand for high- and 
low-index layers). Both effects are totally contrary to 
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FIG. 1. (Color online.) (a) Layout of the experimental set¬ 
up; (b) SEM images of the cross-section of ALD-fabricated 
multilayers with layer thicknesses 10 nm and 20 nm (with 
different layer ordering). (c) Schematics of the four fabricated 
samples differing by the layer thicknesses and ordering. 
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the predictions of the EMA, which is independent of in¬ 
dividual layer features. This anomalous EMA breakdown 
is both enlightening and practically promising (e.g., for 
sensing and switching applications), therefore it would be 
of great importance to observe this effect experimentally. 
However, the effect reported in [H only becomes notice¬ 
able in multilayers containing at least a few hundreds 
of 10-nm thick layers with < 1 nm thickness tolerance, 
which is extremely challenging to fabricate. 

In this Letter, we report the experimental observation 
of the EMA breakdown effect by placing subwavelength 
multilayers on a photonic resonator set-up specially de¬ 
signed to enhance reflectance differences arising from the 
EMA breakdown Q. The resonator makes the effect 
measurable for structures containing only a few tens of 
layers significantly relaxing the fabrication requirements. 
The multilayers themselves were fabricated using atomic 
layer deposition (ALD), chosen for its several advantages 
compared to other deposition techniques (e.g. precise 
thickness control, excellent step coverage, and conformal 
deposition Hi). Our measurements confirm that the re¬ 
flection spectra of the subwavelength multilayers become 
sensitive to both layer thickneses and layer ordering in a 
range of incident angles immediately preceding the TIR 
angle, in full agreement with the theoretical predictions. 

The layout of the experimental arrangement is shown 
in Eig. [IJa). The central element of the set-up is a 
subwavelength multilayer film with total thickness 200 
nm, containing alternating low-index (L) layers made of 
alumina (AI2O3) and high-index (H) layers made of ti- 
tania (Ti02). The layer of silicon nitride (Si3N4) be¬ 
neath the multilayer serves as an index matched layer 
(nsiN — \J ^ resonator required to ob¬ 
serve the EMA breakdown. Since the incident angles 
of interest have to be near-TIR, a high-index ambient 
medium is required. Eor this purpose a semi-cylindrical 
prism made of zinc selenide (nznSe > > ^l) was used. 

All the samples were prepared and assembled in a class 
100 cleanroom. The fabrication of the multilayers was 
performed in a hot-wall ALD system (Picosun R 200 ). 
The precursors used for AI2O3 and Ti02 deposition were 
trimethylaluminum A1(CH3)3 and titanium tetrachloride 
TiCU, respectively (both from Sigma-Aldrich). In both 
processes the oxidant source was deionized water. The 
deposition temperature was 120°C in order to prevent 
the crystal anatase phase transition of Ti02 known to 
occur at temperatures above 150 °C [uj which increases 
the films roughness. The growth rates of AI2O3 and Ti02 
films were determined to be 0.047 and 0.089 nm/cycle re¬ 
spectively (in agreement with previously reported data 


m 


using varying-time deposition with ellipsometric 
characterization of the films thicknesses and refractive 
indices (VASE, J.A. Woollam Co.). 

The index matched resonator layer was fabricated by 
depositing 1040 nm of Si3N4 on 100 mm Si(lOO) wafers 
using low-pressure chemical vapor deposition. The pro¬ 


cess was carried out at 780 ° C with ammonia (NH3) and 
dichlorsilane (SiH2Cl2) as reactive gases. The thickness 
and refractive index of the deposited Si3N4 film were 
tested using spectroscopic ellipsometry. The film was 
carefully inspected for cracks, particles and other defects 
using dark field optical microscopy. The wafer with the 
best-quality Si3N4 coating was selected and cleaved in 
pieces, which were used as substrates for the subsequent 
deposition of AI2O3 /Ti02 multilayers. Before inserting 
each substrate into the ALD reactor, it was placed on 
a Si carrier wafer. Therefore the AI2O3 /Ti02 multi¬ 
layers were grown not only on the Si3N4 layer but also 
on the dummy carrier wafer. After the ALD process, the 
dummy wafer was cleaved, and its cross-section was char¬ 
acterized using scanning electron microscopy (SEM). The 
SEM images reveal high-quality homogeneous, conformal 
coatings, as seen in Eig. [IJb). 

We fabricated four different configurations of multilay¬ 
ers arrangement, shown schematically in Eig. [U^c). Two 
of them comprise 20 alternating Ti02 and AI2O3 lay¬ 
ers with d = 10 nm thickness and different layer order¬ 
ing, i.e., whether the layer closest to the substrate is a 
Ti02 or an AI2O3 layer. The other two samples com¬ 
prise 10 alternating layers with double thickness d = 20 
nm, likewise with different layer ordering. Note that in 
the symbolic representation, where H and L refer to 10 
nm thick titania and alumina layers, respectively, these 
four samples can be unambiguously denoted as [HLHL] ^, 
[LHLH]^ [HHLL]^ and [LLHH]^ [see Eig.^c)]. 

If the EMA was valid, all four samples would be 
homogenized to the same D = 200 nm thick slab of 
anisotropic material with permittivity tensor components 

i=i i=i 
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FIG. 2. (Color online.) Calculated reflectance differences (a) 
(b) AR^J^\X,0) for TE-polarized light at near- 
TIR angles of incidence; (c-d) Same as (a-b) but for TM- 
polarized light. The dashed line shows the spectral position 
of the TIR angle ^tir(A) in presence of material dispersion. 
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in the directions parallel and perpendicular to the layers, 
respectively. Here, N is the number of layers and rij is 
the refractive index of the layer 0 . Therefore, the 
EMA predicts that all samples feature identical reflection 
spectra. Thus, the difference in the reflectances between 
the samples signifies the breakdown of the EMA. 

Quantitatively, this breakdown can be analyzed by de¬ 
termining the following two quantities: 

= -^|<i=10 nm R\d=20 nm? . . 

ARo = — ^L..o 

which reflect the sensitivity of the reflectance spectra to¬ 
wards layer thickness d and layers ordering, respectively. 
Specifically for the fabricated samples, we can introduce 

^^d ^ ~ ^[HLHL]5 — ^[HHLL]5, 

^^^d^ ~ ^[LHLH]5 — ^[LLHHjs; 

= i?[HHLL]5 - ^[LLHH]5, 

= i?[HLHL]5 — ^[LHLH]5. 

Eigure [2] presents the calculated theoretical dependen¬ 
cies of ARd and ARq on the wavelength and angle of 
plane wave incidence for both polarizations of light (TE 
and TM). The calculations were performed by the stan¬ 
dard transfer matrix approach as used in @ with refrac¬ 
tive indices of the materials either measured by spectro¬ 
scopic ellipsometry (AI 2 O 3 , Ti 02 , and Si 3 N 4 ), or ob¬ 
tained from literature fsilico n Il4] . ll5| ) and manufacturer- 
provided data (zinc selenide [l 6 |). Thus, dispersion in all 
materials is accurately taken into account. 

In accordance with our earlier theoretical predictions 
for plane wave incidence Eig.[2]shows a series of char¬ 
acteristic peaks in both ARd and ARq with the shape de¬ 
pendent on polarization. These EMA breakdown peaks 
occur at incidence angles 0 below ^tir(A) (see the dashed 
lines in Eig. m. The breakdown effect is generally more 
pronounced in the TE polarization and is more sensitive 
to the ordering of layers than to the variation of layers 
thicknesses. The values of AR are seen to decrease for 
larger wavelengths. However, the typical peak values of 
\AR\ for the target range A = 610-1610nm are between 
0.05 and 0.25, and can even reach 0.5, which is favorable 
for the measurements. Such high values of what would 
otherwise be a very weak effect are brought about by the 
resonator layer beneath the subwavelength structure 

To observe the EMA breakdown effect experimentally, 
we employed a modified Otto-Kretchmann configuration 
as illustrated in Eig. [TJa). A multilayer sample is placed 
in close proximity to a semi-cylindrical ZnSe prism, which 
is sufficiently high-index to achieve TIR at ^ ^ 50°. The 
light source was a super-continuum broadband laser (Su- 
perK, NKT Photonics A/S, A = 600-2500nm). Its col¬ 
limated output beam was polarized by a double Glan- 
Thompson polariser and focused at the ZnSe-sample in¬ 
terface, using a set of parabolic mirrors. The reflected 




FIG. 3. (Color online.) (a) Schematics of the experimental 
set-up for reflectance measurements. Photos of (b) 20nm tita- 
nia terminated sample [LLHH]^ and (c) lOnm alumina termi¬ 
nated sample [HLHL]^ affixed to the ZnSe prism taken from 
the prism side and showing the appearance of the Newton 
rings, as well as the input and output light beams. 


beam was collected to a multimode fiber using another 
parabolic mirror and led into an optical spectrum anal¬ 
yser (OSA, Yokogawa Electric Gorp.) with the measuring 
range A = 350-1750 nm. 

The sample was attached to the prism with a custom- 
made holder tightened by a small-diameter screw. To 
minimize the air gap, which would have a dramatic influ¬ 
ence on reflectance measurements according to the mod¬ 
elling results (especially for TM polarization), and to re¬ 
duce the risk of dust trapping between the prism and the 
sample, the attachment of the sample was performed in 
the cleanroom. The quality of the attachment was mon¬ 
itored visually controlling the appearance of the Newton 
rings around the location of the screw on the holder as 
the screw was tightened. The optical set-up was aligned 
so that the incident beam was aimed at the innermost 
region of the ring pattern as shown in Eig. 3(b-c) 

Incident angle 0 was varied by a rotating stage with 
precision 0.17°; another rotating stage was used to hold 
the collecting block (mirrors and fiber) for the reflected 
beam [Eig. [3l[a)]. The focused incident beam was found 
to span a range 0 of incident angles 0. It was wavelength- 
dependent and estimated to vary from (j) = 1.3° at 
A = 600 nm to 0 = 3.9° at A = 2 jam. Eor each an¬ 
gle intensity spectra were obtained through averaging 
over 9 independent measurements. The measured spec¬ 
tra were normalized by a reference spectrum obtained by 
averaging 3 above-TIR spectra {0 = 53°, 54°, 55°), since 
R = 1 for 0 > 6>tir- The estimated wavelength-averaged 
error in the reflectance spectra was 4%. 
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FIG. 4. (Color online.) Comparison between experimentally 
measured (dashed lines) and theoretically calculated (solid 
lines) reflectance differences: for (a) TE and (b) TM 

polarization; (c) and (d) AR^J"^ for the TE polariza¬ 

tion. The curves for different incident angles (marked on the 
right) are shifted vertically by 0.4 (TE) and 0.2 (TM) for the 
sake of convenience in comparison. 


In order to compare with the experimental measure¬ 
ments, the calculated theoretical spectra were averaged 
over the wavelength-dependent angle interval 0 corre¬ 
sponding to the beam divergence, as described above. 
Additionally, a mismatch in 0 equal to AO = 0.2° (TE) 
and AO = 0.3° (TM) was introduced to compensate for 
set-up misalignment. These values of AO was determined 
by minimizing the error between measured and exper¬ 
imental spectra in multivariate optimization (see Sup¬ 
plementary Information). The optimization additionally 
confirmed that there is no air gap between the prism and 
the thin-film sample, so an optical contact between them 
was achieved without the use of immersion liquid (im¬ 
possible at such high refractive indices). 

Both experimental and theoretical results for ARq and 
ARd in a variety of cases are shown in Fig. |4]for a range of 
0. We see that AR^^^ in the TE polarization [Fig. lU^a)] 
exhibits small-amplitude ripples across the spectrum for 
lower values of 0 (20-30°), in line with the expectation 


of Fig.[2][a). As 0 approaches ^tir, the EMA breakdown 
becomes stronger, and AR reaches values around 0.4- 
0.5. Above 6>tir, Ai? vanishes since light undergoes TIR 
{R = 1) for both samples. Overall, the measured AR^^^ 
behaves very close to the theoretical predictions. 

For the corresponding AR^^^ in the TM polarization 
[see Fig. lU^b)], the reflectance differences have a lower 
amplitude and fewer characteristic spectral features than 
for the TE case. However, non-zero AR both in the¬ 
ory and in experiment is still apparent. Similar behavior 
as for AR^^\ albeit with proportionally smaller ampli¬ 
tudes, is observed for [Fig. lU^c)], as well as for 

AR^^^ [Fig. |4l[d)]. There is still a good agreement be¬ 
tween theory and experiment. The measured reflectance 
difference exceeds the experimental error (4% average 
and 3-6% depending on the wavelength) and reproduces 
most of the theoretically predicted spectral features. 

In summary, we have experimentally demonstrated the 
effect of the EMA breakdown for all-dielectric multilay¬ 
ers with deeply subwavelength layer thicknesses (1/30 to 
1/160 of the incident light wavelength). We have shown 
that the effect is present even for relatively thin struc¬ 
tures, with the total thickness of 200 nm, also smaller 
than the wavelength. The EMA breakdown manifests 
as the difference in the reflectance spectra of structures 
with different layer thickness (20 nm vs. 10 nm) as well 
as different layer ordering (whether the layer closer to 
the substrate is a high- or a low-index layer), as shown 
in Fig. [He). The measured reflectance difference spec¬ 
tra, reaching values of around 0.5, are in good agreement 
with theoretical transfer matrix calculations (Fig. |4]). 

Our results can be used in ellipsometry of multilayer 
structures, both to correct the existing ellipsometry mod¬ 
els that rely on the EMA, and to devise new models 
specifically based on the measurement of the features re¬ 
lated to the EMA breakdown. 

As regards the applicability of the obtained results to 
sensing (using the high sensitivity of the EMA break¬ 
down to the incident angle and the refractive index be¬ 
hind the multilayer i) , we see that the spectral features 
in the reflectance spectra are not as sharp as theoreti¬ 
cally predicted. This is primarily due to the use of the 
hemicylindrical prism together with a focused beam from 
a broadband light source, causing the light waves that 
reach the sample to have a finite range (p of incident an¬ 
gles 0 and therefore smearing the spectral features. It 
is anticipated that adjusting the experimental set-up to 
reduce this range of angles (e.g. by using a triangular 
prism or a light source with a collimated output) would 
lead to sharper spectral features in the reflectance dif¬ 
ference spectra (see Fig. El suitable for high-sensitivity 
angle and refractive index measurements. 
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